Measurement of the solar neutrino capture rate with gallium metal. 
Ill: Results for the 2002-2007 data-taking period. 
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The Russian-American experiment SAGE began to measure the solar neutrino capture rate with a target 
of gallium metal in December 1989. Measurements have continued with only a few brief interruptions since 
that time. In this article we present the experimental improvements in SAGE since its last published data 
summary in December 2001. Assuming the solar neutrino production rate was constant during the period of 
data collection, combined analysis of 168 extractions through December 2007 gives a capture rate of solar 
neutrinos with energy more than 233 keV of 65.4^3^ (stat) (syst) SNU. The weighted average of the results 
of all three Ga solar neutrino experiments, SAGE, Gallex, and GNO, is now 66. 1 ± 3 . 1 SNU, where statistical and 
systematic uncertainties have been combined in quadrature. During the recent period of data collection a new 
test of SAGE was made with a reactor-produced ^'Ar neutrino source. The ratio of observed to calculated rates 
in this experiment, combined with the measured rates in the three prior ^'Cr neutrino-source experiments with 
Ga, is 0.87 ± 0.05. A probable explanation for this low result is that the cross section for neutrino capture by the 
two lowest-lying excited states in "Ge has been overestimated. If we assume these cross sections are zero, then 
the standard solar model including neutrino oscillations predicts a total capture rate in Ga in the range of 63 SNU 
to 66 SNU with an uncertainty of about 4%, in good agreement with experiment. We derive the current value of 
the neutrino flux produced in the Sun by the proton-proton fusion reaction to be i/i® = (6.0 ± 0.8) x 10"'/(cm^ s), 
which agrees well with the pp flux predicted by the standard solar model. Finally, we make several tests and 
show that the data are consistent with the assumption that the solar neutrino production rate is constant in time. 

PACS numbers: 26.65.+t, 96.6().-j, 95.85.Ry, 13.15.+g 



I. INTRODUCTION 



The SAGE experiment was built to measure the capture rate 
of solar neutrinos by the reaction ^'Ga H- — > ^'Ge -H and 
thus to provide information to aid in understanding the deficit 
of neutrinos observed in the ^^Cl experiment [11, in which 
only about one-third of the solar neutrino capture rate pre- 
dicted by the standard solar model was detected. The feature 
that distinguishes the Ga experiment from all other past or 
present solar neutrino detectors is its sensitivity to the proton- 
proton fusion reaction, p + p ^ d + + Vg, which generates 
most of the Sun's energy. Ga experiments have provided the 
only direct measurement of the current rate of this reaction. 

A full description of the SAGE experiment and the results 
of each measurement from its inception to December 1997 
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was presented in Ref fT]. Part II of this series, although not 
called by this name, described the changes to the experiment 
and gave the results for the period January 1998 to December 
2001 |3|. In Sees, [n] and III of the present article we do the 
same for the six-year period January 2002 to December 2007. 
We then discuss the four neutrino source experiments with Ga 
in Sec. |IV| give the present interpretation of the SAGE results 
in Sec. jvf derive the contemporary value of the neutrino flux 



produced by the proton-proton fusion reaction in Sec. VI and 
present a brief consideration of the question of possible time 
variation in the data in Sec. IVIII 



In addition to SAGE, there also existed a second Ga solar 
neutrino experiment called Gallex. It contained 30 tons of 
gallium in a solution of GaCls and measured the solar neutrino 
capture rate from 1991 to 1997. In 1998 this experiment was 
reconstituted under the name of GNO and it took data until 
2003. We give the results of these experiments and combine 
them with the SAGE data in Seclml 
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II. EXPERIMENTAL PROCEDURES 
A. Overview 

The SAGE experiment is in a dedicated deep-underground 
laboratory excavated into the side of Mt. Andyrchi in the 
northern Caucasus mountains of Russia. The rock overburden 
is equivalent to 4700 m of water and the measured muon flux 
at the location ofthe experiment is (3.03±0.10)xlO"^/(cm^ s). 

The mass of gallium used in SAGE at the present time is 
50 tonnes. It is contained in seven chemical reactors which are 
heated to 30 °C so the gallium metal remains molten. A mea- 
surement of the solar neutrino capture rate begins by adding 
to the gallium a stable Ge carrier. The carrier is a Ga-Ge alloy 
with a known Ge content of approximately 350 jjg and is dis- 
tributed equally among all reactors. The reactor contents are 
stirred thoroughly to disperse the Ge throughout the Ga mass. 
After a typical exposure interval of one month, the Ge carrier 
and ^'Ge atoms produced by solar neutrinos and background 
sources are chemically extracted from the Ga. The final step 
of the chemical procedure is the synthesis of germane (GeH4), 
which is used as a proportional counter fill gas with an admix- 
ture of (90-95)% Xe. The total efficiency of extraction is the 
ratio of mass of Ge in the germane to the mass of initial Ge 
carrier and is typically (95 + 3)%. 

B. Extraction of Ge from Ga 

The extraction procedures from 1990 to 1997 are described 
in Ref. |l2|. At the beginning of 1998 some minor modifica- 
tions were made as described in Ref. |4|. 

Beginning with the December 2005 extraction, the carrier 
used to measure the extraction efficiency was isotopically en- 
riched in either ^^Ge or ^^Ge. At the end of each extrac- 
tion a sample was taken from the final extraction solution and 
this sample was analyzed with an inductively-coupled plasma 
mass spectrometer to determine the fractional content of the 
various Ge isotopes. The efficiency of Ge extraction from the 
Ga metal was then calculated using the method outlined in 
Appendix |A] This procedure for determining the extraction 
efficiency has the advantage that it gives a direct measure of 
any Ge that may enter the sample from unknown sources. 

C. Counting of "Ge 

^^Ge decays to ^'Ga by pure electron capture with a half 
life of 11. 4 days. Two peaks are observed in the proportional 
counter-the K peak at 10.4 keV and the L peak at 1 .2 keV. The 
counter containing the GeH4 from the extraction is placed in 
the well of a Nal detector that is within a large passive shield 
and is counted for a typical period of 6 months. To reduce the 
influence of ^^^Rn, the volume inside the shield around the 
counters is purged with boil-off gas from a dewar filled with 
liquid nitrogen. 

A completely redesigned proportional counter f5\ began to 
be used with the extraction of April 2001 and has been used 



for all but two extractions since that time. In contrast to the 
usual counters with a solid cathode, the cathode of the new 
counters is made from vapor-deposited carbon, thus eliminat- 
ing the usual dead volume behind the cathode. The dead vol- 
ume is further reduced and end effects are nearly eliminated 
by curving inwards the regions of the counter where the cath- 
ode ends. The cathode and anode leads are sealed into the 
counter body with Mo ribbon that makes the counter leak free 
and ensures excellent gain stability. The cathode is so thin that 
the counter body is transparent, making it possible to visually 
inspect all the internal counter parts. 

During 2004-2005 an extensive series of measurements of 
the efficiency of these new counters was made. The methods 
of measurement were described in Ref. |2| and counter fill- 
ings of ^''Ge, ^'Ge, and ■'^Ar were used. The measured volume 
efficiency of the new counters was 96% with a spread in effi- 
ciency of only ±1% for all counters of this type. This should 
be compared with an average volume efficiency of 89% for 
our original counter design. Further, the fraction of events that 
is degraded in energy was found to be significantly less than 
in the old design. These decreases in degraded fraction com- 
bined with the increase in volume efficiency lead to a quite 
dramatic increase in efficiency for these new counters com- 
pared to the old type, approximately 25% in the K peak and 
10% in the L peak. 

Another innovation in the new counter design is that the 
Suprasil counter body is etched in hydrofluoric acid to a thick- 
ness of ~0.2 mm. This permits calibration of the counter with 
our standard ^^Fe source over nearly its entire volume. As 
an undesired side effect, however, the thin body, combined 
with the very thin cathode, makes these counters sensitive to 
low-energy x rays from local radioactivity. To eliminate this 
response, a graded shield consisting of an outer layer of 1 mm 
of Cu and an inner layer of 3 mm of low-background acrylic 
(to absorb Cu x rays) is placed over the counter body during 
measurement with ^'GeH4. 

The pulses from the proportional counter are sent to a fast 
transient analyzer where they are digitized for 800 ns after 
pulse onset at two different gains, one chosen for the L peak 
and the other for the K peak. The transient digitizer serves 
to differentiate fast-rising ^'Ge pulses from generally slower- 
rising background pulses. This can be seen by comparing 
the upper and lower panels of Fig. [T] which show the pulses 
from the 77 extractions that have been measured in the new 
proportional counters. The upper panel is for all events that 
pass the time cuts for Rn (see Sec. |II D| i, are not high- voltage 
breakdown, do not have a Nal coincidence, and occur during 
the first 30 days of counting. The total live time is 1999.8 
days and there are 2063 events. The lower panel of this fig- 
ure shows the 1545 events that occurred between days 100.0- 
130.1 (the same live time duration as in the upper panel). The 
fast-rising ^'Ge events in the L and K peaks are evident in the 
upper panel but missing in the lower panel because the ^'Ge 
has decayed away. 

Aside from replacing some modules that failed, no changes 
were made to the counting system electronics since their de- 
scription in Ref. L2J. 



SAGE results. Part III 



3 




FIG. 1 : Upper panel: Count rate vs energy and rise time for events 
during the first 30 days of counting. Regions where the L and K 
peaks are predicted to occur based on ^'Fe calibrations are shown 
darkened. There are 427 counts in the L-peak region and 287 counts 
in the A'-peak region. The counts in both regions are a combination 
of events from " Ge decay and background. Lower panel: Equivalent 
graph for all events that occurred during an equal live time interval 
beginning at day 100 after extraction. There are 226 counts in the 
L-peak region and 94 counts in the A'-peak region. 



D. Data analysis 

Based on criteria described in fSJ, a group of events is se- 
lected from each extraction that are candidate ^'Ge decays. 
These events are fit to a maximum likelihood function jS), 
assuming that they originate from an unknown but constant- 
rate background and the exponentially-decaying rate of ^'Ge. 
Because only a few ^'Ge counts are detected from each ex- 
traction, a single run result has a large statistical uncertainty 
and thus little significance. 

Several minor changes in the methods of analysis have oc- 
curred since Part II of this series. These include the following: 

• As discussed in Ref. a small fraction of the decays 
of ^^^Rn are occasionally misidentified as pulses from 
^'Ge. To reduce this effect for ^^^Rn located external 
to the counter we now delete all data that is acquired 
within 2.6 hours after counting begins. In our initial 
analysis this time cut was for only one hour. 



To reduce the influence of ^^^Rn that may enter the 
counter when it is filled, we delete all data from 15 min 
before an event that saturates the energy scale to 3 h 
after each saturated event. The SAGE measurements 
before September 1992, however, were measured in 
counting systems that did not have the capability to rec- 
ognize saturated events |2|. To reduce the number of 
these false ^'Ge events produced by ^^^Rn in these early 
runs we determine for all subsequent runs the differ- 
ence in capture rate in the K peak between the data an- 
alyzed with and without this time cut and then subtract 
this difference from the result of each of the runs before 
September 1992. 

The predicted position and resolution of the L peak 
from calibrations was changed slightly from previous 
work. These changes were indicated by the results of a 
set of new calibrations made with counters filled with 
^'Ge. 

The L- and K-peak shapes were changed from pure 
Gaussian to Gaussian plus a degraded term. The new 
functional form for the line shape as a function of en- 
ergy is 



F(E) = he 



-[(£-C)/(V2o-)]^ 



+ hd 



n cr 
2C 



erfc 



V2cr^ 



(1) 



where h, C, and cr are the peak height, center, and width 
and li is a parameter related to the fraction of degraded 
events. The error function term here is the integral of 
the Gaussian from energy E to oo; it is essentially flat 
below the peak, monotonically decreases in the peak re- 
gion, and is zero above the peak. This new line shape 
only makes a very small change to the counting effi- 
ciency in the L peak for a few runs whose energy win- 
dow width is obliged to be less than 2 full widths at half 
maximum. 

• For all runs after August 1992 the likelihood function 
was modified to include a factor that weights each event 
according to its measured energy. This requires knowl- 
edge of the energy distribution for ^'Ge pulses and for 
background events, both of which can be determined 
from the long duration of counting data that we have 
accumulated. When this method is applied, it is found 
that the overall statistical uncertainty decreases by (0.1- 
0.2) SNU, but the systematic uncertainty increases by 
-0.1 SNU. 

These changes in analysis methods have been applied to all 
data. 



III. RESULTS 

^'Ge has been extracted from the Ga target to measure the 
solar neutrino capture rate every month from January 2000 to 
the present time. We even were able to make six solar extrac- 
tions during the time of the ^^Ar neutrino source experiment 
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FIG. 2: Combined SAGE results for each year Shaded band is the combined best fit and its uncertainty for all years. Vertical error bars are 
statistical with 68% confidence. 



TABLE I: Summary of known systematic effects and their uncertain- 
ties. SNU values for extraction and counting efficiency are based on 
a rate of 65.4 SNU. 



Uncertainty 



Origin of uncertainty 


in percent 


in SNU 


Extraction efficiency 






Ge carrier mass 


±2.1% 


±1.4 


Mass of extracted Ge 


+2.5% 


±1.6 


Residual Ge carrier 


±0.8% 


±0.5 


Ga mass 


±0.3% 


±0.2 


Total (extraction) 


±3.4% 


±2.2 


Counting efficiency 






Volume efficiency 


±1.0% 


±0.7 


End losses 


±0.5% 


±0.3 


Monte Carlo interpolation 


±0.3% 


±0.2 


Shifts of gain 


-1.1% 


+0.7 


Resolution 


+0.5%, -0.7% 


-0.3, +0.5 


Rise time limits 


±1.0% 


±0.7 


Lead and exposure times 


±0.8% 


±0.5 


Total (counting) 


+ 1.8%, -2.1% 


-1.2,+1.4 


Nonsolar neutrino production of 


7'Ge 




Fast neutrons 




< -0.02 






< -0.04 


22«Ra 




< -0.7 


Cosmic-ray muons 




< -0.7 


Total (nonsolar) 




< -1.0 


Background events that mimic " 


Ge 




Internal ^-^Rn 




< -0.2 


External ^^^Rn 




0.0 


Internal "^'Ge 




< -0.6 


Total (background events) 




< -0.6 


Energy weighting in analysis 




±0.1 


Total 




-2.8, +2.6 



in 2004 by sending the samples to Gran Sasso. In a coopera- 
tive effort f?!, the GNO collaboration synthesized GeH4 and 
measured the samples in their counting system. 

The results for each individual extraction are tabulated in 



Appendix |B] and the combined result of each year of SAGE 
data since its beginning is shown in Fig.|2] 

The systematic uncertainties in the experiment have been 
considered in detail in Ref. 12] |3] and the most recent val- 
ues are given in Table |l] The only significant changes from 
our previous articles are due to the new proportional counters. 
Their high stability and efficiency have led to a considerable 
reduction of the uncertainties associated with counting. 

In radiochemical experiments the capture rate has been con- 
ventionally expressed in "SNU units", defined as one neu- 
trino capture per second in a target that contains 10"''' atoms 
of the neutrino-absorbing isotope, in our case ^'Ga. For 
all SAGE data from January 1990 through December 2007 
(168 runs and 310 separate counting sets) the global best fit 
capture rate is 65.4;|;2 o SNU, where the uncertainty is statis- 
tical only. If one considers the L-peak and /T-peak data sep- 
arately, the results are 67.2+4-^ SNU and 64.0+^(1 SNU, re- 
spectively. The agreement between the two peaks serves as 
a strong check on the robustness of the event selection crite- 
ria. Including the systematic uncertainty, our overall result is 
65.4+i^;i (stat)!^-6 (syst) SNU. 

As further evidence that we are truly counting ^'Ge, we can 
allow the decay constant during counting to be a free variable 
in the maximum likelihood fit, along with the combined ^'Ge 
production rate and all the background rates. The best fit half- 
life to all selected events in both L and K peaks is then 11.5 + 
0.9 (stat) days, in agreement with the measured value |8J of 
11.43 + 0.03 days. 

The waveform data from the Gallex experiment has re- 
cently been re-evaluated by Kaether using a new pulse-shape 
analysis method ID and the result is 73.1+^,'+^^ SNU. The re- 
sult of the GNO experiment was 62.9+^-^+^'^ SNU | la]. If we 
combine the statistical and systematic uncertainties in quadra- 
ture, then the weighted combination of all the Ga experiments, 
SAGE, Gallex, and GNO, is 

66.1 ± 3.1 SNU. (Present Ga experiment result.) (2) 
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FIG. 3: Results of all neutrino source experiments with Ga. Gallex 
results are from the recent pulse shape analysis of Kaether 1 9 1 ; SAGE 
results are from Refs. 1 1 1 1 and |4J. Hashed region is the weighted 
average of the four experiments. 



IV. SOURCE EXPERIMENTS 

The experimental procedures of the SAGE and Gallex ex- 
periments, including the chemical extraction, counting, and 
analysis techniques, have been checked by exposing the gal- 
lium target to reactor-produced neutrino sources whose ac- 
tivity was close to 1 MCi. SAGE has irradiated about 25% 
of their target with a ^'Cr source LI IJ and an ^^Ar source 
||4][T2][T3] and Gallex has twice used ^'Cr sources to irradiate 
their entire target [ 14 1. The results, expressed as the ratio R of 
the measured ^'Ge production rate to that expected due to the 
source strength, are shown in Fig. |3] The weighted average 
value of the ratio for the four experiments is R - 0.87 + 0.05, 
more than two standard deviations less than unity. Although 
the distribution of results is somewhat unusual, with none of 
the central values from the four measurements lying within 
the Icr band around the weighted average, the quality of fit 
to the average value is quite high (y^/DOF = 1.9/3, GOF - 
59%). 

We can suggest several possibilities for the unexpectedly 
low result in the source experiments: 

1 . We do not correctly know the various efficiency factors 
that enter into the calculation of the production rate, 
namely the extraction efficiency and the counting effi- 
ciency. Both SAGE and Gallex have, however, made 
many ancillary experiments ||3] [T4ll that have estab- 
lished with high probability that these efficiencies and 
their accompanying systematic uncertainties are well 
determined. These tests have also proven that there are 
no substantial errors in the methods used to select ^'Ge 
events or in the methods of analysis. Fuither, the ^'As 
experiment of Gallex [15 1 has ruled out any "hot-atom" 
chemical effects that might make the ^'Ge atoms pro- 
duced by neutrino capture difficult to extract. We thus 
very strongly doubt that the low average result of the 
source experiments is due to incorrect knowledge of ef- 
ficiencies, errors in event selection, improper function- 
ing of ffie counting systems, or errors in analysis. 



2. A statistical fluctuation. A test of the compatabil- 
ity of the four source experiments to R = 1.0 gives 
X^/GOF - 7.7/3, whose probability is 5.3%. The prob- 
ability is small, but still quite possible. 

3. Electron neutiinos disappear due to a real physical ef- 
fect of unknown origin. Some possibilities that have 
been suggested are a transition to sterile neutrinos lfT6ll 
or quantum decoherence in neutrino oscillations IfTTl . 

4. The production rate from the source is not as great as 
has been assumed. It is our opinion that this is ffie 
most likely cause of the apparently low result in the 
source experiments. As suggested by Haxton |18 | it 
is quite possible that the cross sections for neutrino 
capture to the two lowest excited states in ^'Ge, both 
of which can be reached using either ^'Cr or ^^Ar 
sources, have been overestimated. 95% of the capture 
rate with these sources arises from the ^'Ga to ^'Ge 
ground-state transition with 5% due to transitions to 
the two excited states. If the contribution of the ex- 
cited states to the predicted rate were to be zero then 
R = /?(measured)//?(predicted) - 0.92 ± 0.06 and the fit 
to the expected value of 1 .0 becomes quite reasonable 
(X^/DOF = 4.58/3, GOF = 21%). 

A concern in this context is the implication of the appar- 
ently low result of the source experiments on the solar neu- 
trino result given in Eq. It is difficult to address this con- 
cern because we do not understand why the source experi- 
ments give a lower result than expected. If we suppose that the 
cause is item 1 in the list above, then the rate in Eq. Q should 
be divided by the factor 0.87, i.e., we should add 15% to the 
systematic uncertainty. But, as stated above, we consider that 
explanation for the apparent discrepancy in the source exper- 
iments to be very unlikely. However, if we suppose that the 
cause of the low result in the source experiments is any of 
the other items in the list above, then the source experiments 
have no bearing on the solar neutrino result and the rate in 
Eq. Q should not be changed. Because we do not know why 
the source experiments appear to be low, we can only caution 
the reader to accept the result in Eq. (|2]i on a provisional ba- 
sis, subject to the caveats that not all effects in the emission of 
neutrinos from the Sun and ffie capture of neutrinos by ^'Ga 
may be fully understood. 

Based on the information given in the definitive article of 
Bahcall ||T9l on ffie neutrino capture cross section of ^'Ga, we 
have approximately calculated the cross section if we assume 
zero strength for capture to the first two excited states of ^'Ge. 
These results are given in Appendix |C] and will be used as a 
working hypothesis in what follows. 



V. INTERPRETATION OF RESULTS 

In contrast to all other present or past solar neutrino exper- 
iments, the radiochemical Ga experiment, because of its low 
threshold of 233 keV, is sensitive to all components of ffie 
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TABLE II: Solar neutrino fluxes calculated by the standard solar 
model of Bahcall, Pena-Garay, and Serenelli 1211 for two dilferent 
conservative choices of heavy element composition, labeled GS98 
1 22 1 (high metallicity) and AGS05 |23| (low metallicity). The spec- 
trum components refer to the nuclear reaction from which they origi- 
nate. Units of flux are 10'" (pp), 10" ('Be), 10** (pep, "N, ''O), 10* 
(^B, '^F), and 10^^ (hep) cm^^s"'. The uncertainty values are at 68% 
confidence. 



Spectrum 




Flux (pf 




component i 


GS98 




AGS05 



PP 


5 97(1+0.006-, 


5 04(1 +0.005-, 


pep 


1 41(1+"'"') 


i-45(it!!l!J) 


'Be 


5.07(1!»:«^) 


4-55(1^!!^) 


13N 


2.88(1!0|5) 


i.89(i!»:W) 




2-15(1^1^) 


i.34(i!°:[«) 


Hp 


5.82(1!«:«) 


3.25(1!°:|«) 


**B 


5.94(rO;]|) 


4.72(1+0;]|) 


hep 


7.90(1!»:|^) 


8.22(i:»;[5) 



component as the neutrinos are produced at different locations 
in the Sun and thus pass through regions of different electron 
density during their travel to the Earth. Pf{E) can only be cal- 
culated if one knows where in the Sun the neutrinos are made 
and thus requires a solar model. 

We integrate Eq. (jSjl and obtain A, = where 



{PT) 



f 

Jo 



Sf{E)Pl\E)dE 



(6) 



is the spectrum-weighted average value of P^*". The physical 
interpretation of {Pf) is as the ratio of the solar neutrino am- 
plitudes at the surface of the Earth and at the production point 
in the Sun: 



{PT) 



A©' 



(7) 



Combining these equations, and expressing R as the sum of 
its spectral components, R - Xi, /?,, we have 



solar spectrum, from the low-energy pp neutrinos to the high- 
energy neutrinos produced in the decay of ^B. In Table [ll] we 
give the flux of the various solar neutrino components at their 
production regions in the Sun as calculated by Bahcall and 
collaborators |20l |2T|. In this section we will estimate the 
neutrino capture rate from each flux component and compare 
their total to the measured rate. 

The total capture rate R of solar neutrinos in a radiochemi- 
cal experiment such as Ga is given by 



R 



I cr(£)(D* (£)£/£, 

•^£'llir,-^h,il,l 



(3) 



where cr(E) is the cross section of the neutrino-capture reac- 
tion and is the total flux of electron neutrinos at the 
Earth, which can be expressed as 



(4) 



In this expression the index / refers to the various nuclear re- 
actions in the Sun that produce neutrinos [pp, 'Be, pep, ^^N, 

'^O, ''F, ^B, and hep), (f>f is the amplitude of flux component 

/ at the Earth, and sf(E) is the spectrum of the /th neutrino 
component at the Earth, each of which is normalized such that 
£° sf{E)dE - 1. The neutrino spectrum at the Earth is re- 



lated to the spectrum produced in the Sun S f'iE) by 
sf{E)^AiSf{E)P^'{E), 



(5) 



where A,- is a constant of normalization and P^iE) is the prob- 
ability that an electron neutrino produced in the Sun by reac- 
tion / with energy E will reach the Earth without a change of 
flavor, commonly called the survival factor. The physical ori- 
gin for the reduction of the electron component of the solar 
neutrino flux is the now well-established mechanism of MS W 
neutrino oscillations l,24J . P^iE) is different for each flux 



where 



/-*oc 

I o-{E)Sf{E)P^\E)dE, 



(8) 



(9) 



or, equivalently, if it is the flux at the Earth that is assumed 
known. 



where 



Ri 



(10) 



(11) 



In Table 



III 



values of {P^) and {erf) for 
These were calculated assum- 



we give 

each neutrino~component 
ing three-neutrino mixing to active neutrinos with parame 



ters from Ref. 



Am?, 



(7.65+;}23) X 



0.20^ 



33.46+! nfj degrees, and - 5.7+^-^ degrees 



lO-=eV^0l2 = 
The approxi- 



-i.oo 

mate formulae given in Ref. Il26l were used for the survival 
probability PfiE). As we show in Sec. 



VII 



there is no ap- 
preciable difference between the day and night capture rates 
in Ga and thus regeneration in the Earth was neglected. The 
cross sections cr{E) were taken from Appendix [C] for Ga 
and Ref. L27J for CI. The neutrino spectra <pf(_E) are from 
Refs. m {pp, ^^N, '^O, "F), E?) ('^B), and 128) {hep). 

Now that all the terms have been calculated, we can use 
the fluxes in Table |Ilj combined with Eqs. ([8]l and (111, to 
predict the capture rate in Ga from each of the solar neutrino 
components. The individual rates and the total rate are given 
in Table |IV] for two recent solar models from Table |ll] For 
both models there is good agreement between the calculated 
total rate and the observed capture rate of 66.1 + 3.1 SNU. 
The major contribution to the uncertainty in the predicted total 
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TABLE ni: Factors needed to compute the solar neutrino capture rate in ^'Ga and ''^Cl solar neutrino experiments. The uncertainty values are 



at 68% confidence. The parameter (crp) is defined in Eq. 111'. 





S ect 
pec . 




Percent uncertainty in (P": 


due to 


Total unc. 


W i 1 




Percent uncertainty 


in {a^ 


ue o 






1 

loi^i Line. 


Exp. 


comp. 






f 12 


"13 




in (Pf > (%) 


(lO"-"" cm-) 


cr 






012 




fil3 




in (a-f) (%) 




PP 


0.561 


+ O.O.- 0.0 


+ 2.2.- 2.8 


+ 2.0,- 


3.1 


+ 3.0,- 


4.2 


11.75 


+ 2.4,- 


2.3 


+ 0.0,- 0.0 


+ 0.0,- 


0.0 


+ 0.0,- 


0.0 


+ 2.4,- 2.3 




pep 


0.521 


+ 0.3,- 0.3 


+ 1.9,- 2.3 


+ 1.9,- 


2.9 


+ 2.7,- 


3.7 


194.4 


+ 17 ,- 


2.4 


+ 0.0,- 0.0 


+ 0.0,- 


0.0 


+ 0.0,- 


0.0 


+ 17 ,- 2.4 




'Be 


0.542 


+ 0.1,- 0.1 


+ 2.1,- 2.6 


+ 2.0,- 


3.0 


+ 2.9,- 


4.0 


68.21 


+ 7.0,- 


2.3 


+ 0.0,- 0.0 


+ 0.0,- 


0.0 


+ 0.0,- 


0.0 


+ 7.0,- 2.3 




"N 


0.545 


+ 0.1,- 0.1 


+ 2.1,- 2.6 


+ 2.0,- 


3.0 


+ 2.9,- 


4.0 


56.83 


+ 9.8,- 


2.3 


+ 0.0,- 0.0 


+ 0.1,- 


0.0 


+ 0.0,- 


0.0 


+ 9.8,- 2.3 




'50 


0.535 


+ 0.2,- 0.2 


+ 2.0,- 2.5 


+ 1.9,- 


3.0 


+ 2.8,- 


3.9 


107.1 


+ 13 ,- 


2.3 


+ O.O.- 0.0 


+ 0.1.- 


0.1 


+ 0.0,- 


0.0 


+ 13,- 2.3 




"F 


0.535 


+ 0.2,- 0.2 


+ 2.0,- 2.5 


+ 1.9,- 


3.0 


+ 2.8,- 


3.9 


107.7 


+ 13,- 


2.3 


+ 0.0,- 0.0 


+ 0.1,- 


0.1 


+ 0.0,- 


0.0 


+ 13,- 2.3 




«B 


0.365 


+ 0.7.- 0.6 


+ 3.2,- 2.2 


+ 1.8,- 


2.7 


+ 3.8,- 


3.5 


21400 


+32, 


14 


+ 0.2,- 0.2 


+ 2.1,- 


1.7 


+ 0.1,- 


0.1 


+32 ,-15 




hep 


0.337 


+ 0.5,- 0.5 


+ 4.8,-3.4 


+ 1.8,- 


2.8 


+ 5.2,- 


4.4 


66000 


+33 , 


15 


+ 0.2,- 0.2 


+ 1.4,- 


1.1 


+ 0.1,- 


0.1 


+33 ,-16 




pep 


0.521 


+ 0.3,- 0.3 


+ 1.9,- 2.3 


+ 1.9,- 


2.9 


+ 2.7,- 


3.7 


16.00 


+ 2.O.- 


2.0 


+ 0.0,- 0.0 


+ 0.0,- 


0.0 


+ 0.0,- 


0.0 


+ 2.0,- 2.0 




'Be 


0.542 


+ 0.1,- 0.1 


+ 2.1,- 2.6 


+ 2.0,- 


3.0 


+ 2.9,- 


4.0 


2.397 


+ 2.0,- 


2.0 


+ 0.0,- 0.0 


+ 0.0,- 


0.0 


+ 0.0,- 


0.0 


+ 2.0,- 2.0 




"N 


0.545 


+ 0.1,- 0.1 


+ 2.1,- 2.6 


+ 2.0,- 


3.0 


+ 2.9,- 


4.0 


1.686 


+ 2.0,- 


2.0 


+ 0.1,- 0.1 


+ 0.1,- 


0.1 


+ 0.0,- 


0.0 


+ 2.0,- 2.0 




'50 


0.535 


+ 0.2,- 0.2 


+ 2.0,- 2.5 


+ 1.9,- 


3.0 


+ 2.8,- 


3.9 


6.662 


+ 2.0,- 


2.0 


+ 0.1,- 0.1 


+ 0.2,- 


0.1 


+ 0.0,- 


0.0 


+ 2.0,- 2.0 




"F 


0.535 


+ 0.2,- 0.2 


+ 2.0,- 2.5 


+ 1.9,- 


3.0 


+ 2.8,- 


3.9 


6.710 


+ 2.0,- 


2.0 


+ 0.1,- 0.1 


+ 0.2,- 


0.1 


+ 0.0,- 


0.0 


+ 2.0,- 2.0 




«B 


0.365 


+ 0.7,- 0.6 


+ 3.2,- 2.2 


+ 1.8,- 


2.7 


+ 3.8,- 


3.5 


I0I40 


+ 3.7,- 


3.7 


+ 0.2,- 0.2 


+ 2.4,- 


1.9 


+ O.I,- 


0.1 


+ 4.4,- 4.1 




hep 


0.337 


+ 0.5,- 0.5 


+ 4.8,-3.4 


+ 1.8,- 


2.8 


+ 5.2,- 


4.4 


40910 


+ 3.7.- 


3.7 


+ 0.2,- 0.2 


+ 1.5.- 


1.2 


+ 0.1,- 


0.1 


+ 4.0,- 3.9 



TABLE IV: Capture rates Rj for Ga experiments calculated with fluxes from Ref. |21 









With GS98 composition 














With AGS05 composition 










Spect. Cap. rate 




Percent uncertainty in rate due to 






Total unc. 


Cap. rate 




Percent uncertainty in rate due to 






Total unc. 


comp. 


(SNU) 


•P 




AmJ, 






ei3 




in rate {%) 


(SNU) 







Am^j 






^13 




in rate (%) 


PP 


39.35 


+ 0.6,- 0.6 


+ 2.4,- 2.3 


+ 0.0,- 0.0 


+ 2.2,- 


2.8 


+ 2.0,- 


3.1 


+ 3.9,-4.8 


39.81 


+ 0.5,- 0.5 


+ 2.4,- 2.3 


+ 0.0,- 0.0 


+ 2.2,- 


2.8 


+ 2.0,- 


3.1 


+ 3.9,- 4.8 


pep 


1.43 


+ 1.1,- I.l 


+ 17.0,- 2.4 


+ 0.3,- 0.3 


+ 1.9,- 


2.3 


+ 1.9,- 


2.9 


+ 17.2,- 4.6 


1.47 


+ 1.0,- 1.0 


+ 17.0,- 2.4 


+ 0.3,- 0.3 


+ 1.9,- 


2.3 


+ 1.9,- 


2.9 


+ 17.2,- 4.5 


'Be 


18.73 


+ 6.0,- 6.0 


+ 7.0,- 2.3 


+ 0.1,- 0.1 


+ 2.1,- 


2.6 


+ 2.0,- 


3.0 


+ 9.7,- 7.5 


16.81 


+ 6.0,- 6.0 


+ 7.0,- 2.3 


+ 0.1,- 0.1 


+ 2.1,- 


2.6 


+ 2.0,- 


3.0 


+ 9.7,- 7.5 


13n 


0.89 


+ 15.0,-15.0 


+ 9.8,- 2.3 


+ 0.1,-0.1 


+ 2.1,- 


2.6 


+ 2.0,- 


3.0 


+18.1,-15.7 


0.58 


+ 14.0,-13.0 


+ 9.8,- 2.3 


+ 0.1,- 0.1 


+ 2.1,- 


2.6 


+ 2.0,- 


3.0 


+ 17.3,-13.8 


"0 


1.23 


+ 17.0,-16.0 


+ 12.9,- 2.3 


+ 0.2,- 0.2 


+ 2.0,- 


2.4 


+ 1.9,- 


3.0 


+21.5,-16.6 


0.77 


+ 16.0,-15.0 


+ 12.9,- 2.3 


+ 0.2,- 0.2 


+ 2.0,- 


2.4 


+ 1.9,- 


3.0 


+20.7,-15.6 


"F 


0.03 


+ 19.0,-17.0 


+ 12.9,- 2.3 


+ 0.2,- 0.2 


+ 2.0,- 


2.4 


+ 1.9,- 


3.0 


+23.1,-17.6 


0.02 


+ 16.0,-15.0 


+ 12.9,- 2.3 


+ 0.2,- 0.2 


+ 2.0,- 


2.4 


+ 1.9,- 


3.0 


+20.8,-15.6 


*B 


4.64 


+ 11.0,-11.0 


+31.8,-14.4 


+ 0.5,- 0.4 


+ 5.4,- 


3.9 


+ 1.8,- 


2.8 


+34.1,-18.7 


3.68 


+ 11.0,-11.0 


+31.8,-14.4 


+ 0.5,- 0.4 


+ 5.4,- 


3.9 


+ 1.8,- 


2.8 


+34.1,-18.7 


hep 


0.02 


+ 15.0,-15.0 


+32.7,-15.4 


+ 0.3,- 0.3 


+ 6.2,- 


4.5 


+ 1.9,- 


2.9 


+36.5,-22.2 


0.02 


+ 15.0,-15.0 


+32.7,-15.4 


+ 0.3,- 0.3 


+ 6.2,- 


4.5 


+ 1.9,- 


2.9 


+36.5,-22.2 


Total 


66.31 


+ 1.9,- 1.9 


+ 3.3,- 1.8 


+ 0.1,-0.1 


+ 1.5,- 


1.8 


+ 1.3,- 


2.0 


+ 4.3,-3.8 


63.16 


+ 1.8,- 1.8 


+ 3.1,- 1.8 


+ 0.1,- 0.0 


+ 1.5,- 


1.9 


+ 1.4,- 


2.1 


+ 4.1,-3.8 



rate is from the neutrino capture cross section, with smaller 
contributions from the solar model flux, 612, and ^n. 

In this analysis we have used the cross sections in Ap- 
pendix [Cj in which the contribution of the two lowest-lying 
excited states in ^'Ge has been set to zero. If instead we use 
the original Bahcall cross sections, then the total rate increases 
by 1.2 SNU with the GS98 composition and by 1.1 SNU with 
the AGS05 composition. Whatever cross sections are assumed 
is thus not a significant factor in the interpretration of the total 
rate in the Ga experiment. The cross sections are, however, of 
vital importance in understanding the origin of the unexpect- 
edly low result in the source experiments. 

The attentive reader may be concerned that there is a log- 
ical inconsistency in the argument presented here: the pre- 
dicted capture rates we derive for the Ga experiment depend 
on the neutrino oscillation parameters, but the measured total 
rate of the Ga experiment is itself one of the inputs used to de- 
termine the oscillation parameters. Although this is true in a 
strict sense, the neutrino oscillation parameters derived from 
a global fit of all experiments are for all practical purposes 
independent of the rate in the Ga experiment. Rather, the pa- 
rameter 012 is principally determined by the SNO experiment, 
013 by the CHOOZ experiment, and Amjj by the KamLAND 
experiment. Although it was not true in the past, the result of 
the Ga experiment is at present only a very minor input to the 
deteiTnination of these parameters. 

Incidentally, if we carry out the same analysis for the 



CI solar neutrino experiment, the total calculated rate 
is 3.09(1!°°^4) using fluxes based on GS98 and 

2.53(1+°°^^) SNU using fluxes based on AGS05. These 
should be compared with the experimental rate of 2.56 + 
0.16 (Stat) + 0.16 (syst) SNU ffl. 

VI. THE PP NEUTRINO FLUX FROM THE SUN 

In this section we will use the Ga measurement given in 
Eq. (|2]l and the results of other solar neutrino experiments to 
determine the pp flux from the Sun. The conventional way to 
make this calculation is by a combined fit to all experiments, 
as for example is presented in Ref. [291 and fW|. Here we 
give an alternate approach that successively decomposes the 
total measured rate into the components from each neutrino 
source. The final result is identical to what one obtains in a 
combined fit and has the advantage that the argument is simple 
and transparent. 

The rate in Eq. (j2| is the sum of the rates from all the com- 
ponents of the solar neutrino flux, which we denote by 

[pp+'^Be+CNO+pep+^BlGa] = 66.1(1 ± 0.047) SNU. (12) 

We ignore the tiny hep contribution and combine the ^^N, '^O, 
and '^F components into a single value, called here "CNO". 

In an experiment of great technical difficulty, the ^Be flux 
has been directly measured by Borexino and they report the 
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result as 



Using Eqs. dSb and ( 1 1 



5.18(1 ± 0.098) X 10'* neutrinos/(cm2 s) f30\. 



we multiply this flux by the electron 
neutrino survival factor for ^Be and by the cross section of 
^Be on Ga (the values of these factors and their uncertainties 



are given in Table 111 i and obtain the rate of ^Be in Ga of 



fBelGa] = 19.1(l!|;:!J) SNU. 



(13) 



The B flux at the Earth has been directly measured by SNO 
to be = ( 1 .67 ± 0.05) x 10*" electron neutrinos/(cm^ s) ([31] . 
In a similar way to ^Be, we multiply this flux by the spectrum- 
integrated cross section for ^B neutrinos on Ga and obtain the 
**B contribution to the Ga experiment of 



[^B|Ga] = 3.6(l!°:J2) SNU. 



(14) 



Subtracting these measured rates of ^Be and **B from the 



total Ga rate in Eq. ( 12 1 gives 



[pp+CNO+pep\Ga] = 43.3(l!^:[|^;j) SNU. 



(15) 



We can obtain an approximate value for the contribution 
of CNO and pep to the Ga experiment from the measured 
capture rate in the CI experiment [^Be+CNO+/?e;?+'*B|Cl] - 
2.56(1 ± 0.088) SNU |1 1. As in the case of Ga, we use the 
^Be flux measured by Borexino, the ^B flux measured by 
SNO, and the cross sections in Table 111 to determine [^Be|Cl] 



= 0.67(1+[; |[;5) SNU and [**B|C1] = 1.73(1![;°^^) SNU. We 
subtract these values from the total CI rate and are left with 
[CNO+pep\Cl] = 0.19(1![;^^) SNU. 

If we attribute this entire rate to the neutrinos from pep 
then, using the cross sections for pep on CI and Ga, we cal- 
culate a rate of [peplGaJtest = 2.35(l+j Q„) SNU. On the other 
hand, if we attribute this entire rate to CNO, we obtain in the 
same manner a rate of [CNO|Ga],e,st = 3.11(l+j-^^) SNU. The 
upper extreme of these two test rates is 3.11 x (1 H- 1.37) - 
7.37 SNU. As a reasonable estimate we can thus set the sum 
of CNO and pep rates at half this upper limit with an uncer- 
tainty of 100%: 



[CNO-i-peplGa] = 3.68(l!!:J^g) SNU. 



(16) 



We subtract this estimate for the CNO plus pep rate from 



the rate in Eq. ( 15 1 and obtain the result for the measured pp 
rate in the Ga experiment 



[pp\Ga] = 39.7(1 snu. 



(17) 



Dividing this capture rate by the cross section for capture of 



pp neutrinos from Table 111 gives the measured electron neu- 
trino pp flux at Earth of 



4 = 3.38(l!°:[^) X 10'"/(cm2 s). 



(18) 



If we use Eq. and the value of <Pf ) = 0.561(1+JJ JJ^°) from 
Table 111 then tne pp flux produced in the Sun is 



6.0(1+ 0.14) xlO'Vcm^s). 



(19) 



Our present result for the pp flux is in good agreement with 
the previous estimates that we have made during the past six 



years |l3l [32ll33]| , with the major change being a gradual reduc- 
tion of the uncertainty. In the future, as Borexino continues to 
collect data, and as direct measurements are made of the CNO 
and pep fluxes, the uncertainty in this flux should be further 
reduced and eventually may be dominated by the uncertainty 
in the Ga rate itself. By that time, however, there will hope- 
fully be direct experiments that measure the pp flux in real 
time. 

For comparison, we see from Table that the predicted pp 
flux from the two recent solar models with different compo- 
sition is cpfp = 5.97 ± 0.04 and 6.04 ± 0.03, both in units of 
lO'" V(./(cm^ s). There is good numerical agreement between 
these flux values and the result in Eq. ( 19 1, but, as made clear 



by Bahcall and Pena-Garay |34|, there is a large difference in 



interpretation: the result in Eq. ( 19 1 was derived from contem- 
porary solar neutrino experiments and is the pp neutrino flux 
at the present time. In contrast, energy generation in the solar 
model is highly constrained by the measured solar luminosity 
and thus, when the luminosity constraint is imposed, as is the 
case for the models given in Table [H] the calculated pp flux 
is what the Sun was producing some 40 000 years ago. The 
agreement between the present pp flux, as measured by the Ga 
experiment, and the past flux, as inferred from the solar model 
with the luminosity constraint, implies that the pp flux from 
the Sun has not altered (within our 14% uncertainty) during 
the last 40 000 years. 



VII. CONSIDERATION OF TIME VARIATION 

In a plot of the SAGE results as a function of time there is 
a slight visual hint of a long-term decrease, as illustrated in 
Fig. |2] The average rate prior to 1996 is somewhat higher 
than after 1996. A plot of the Gallex-GNO data shows a 
similar behavior [lO]. When examined quantitatively, how- 
ever, the evidence for a long-term decrease in the capture rate 
is unconvincing. A test applied to these yearly SAGE 
data points assuming the rate is constant at 65.4 SNU gives 
X^I\iO¥ = 12.0/17, which has a probability of 80%. The fit to 
a constant rate is thus quite good. 

In previous articles we have demonstrated the agreement 
between the assumption of a constant production rate and the 
SAGE measurements by use of the cumulative distribution 
of the capture rate C{p), defined as the fraction of data sets 
whose capture rate is less than p. Figure |4] shows this distri- 
bution for the data and the expected distribution derived from 
100 simulations of all 168 runs, where it is assumed in the sim- 
ulations that the production rate is constant and has a value of 
65.4 SNU. For each run the rates from the separate L and K 
peaks are used in this figure, not the rate from the L+K combi- 
nation. To ensure that the simulations parallel the real data as 
closely as possible, all parameters of the simulation, such as 
background rates, efficiencies, exposure times, and counting 
times, were chosen to be the same as for the real data. Only 
the number of counts in each run and the times when these 
counts occurred were allowed to vary. 

The data spectrum and the simulated spectrum are very 
similar to each other, indicating that the distribution of cap- 
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FIG. 4: Measured capture rate for all 310 SAGE data sets (jagged 
curve) and the expected distribution derived by 100 Monte Carlo sim- 
ulations of each set (smooth curve). 
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FIG. 5: Lomb power spectrum from all 168 SAGE data runs. The 
mean time of exposure was used as the time of measurement. 



ture rates is what one would expect if the rate is constant. A 
quantitative comparison can be made by calculating the Nw^ 
test statistic for the data distribution and comparing it to the 
distribution from simulations using the method described in 
Ref. |35|. The fraction of simulated spectra whose Nw^ was 
larger than for the data distribution is 43%, which shows that 
the assumption of a constant capture rate is in good agreement 
with our measurements. 

A standard method to look for periodic signals in unevenly 
sampled data, such as we have in SAGE, was devised by Lomb 
and Scargle. Application of this method, using the implemen- 
tation of Press et al. |36|, to all runs from the SAGE exper- 
iment yields the power spectrum shown in Fig. |5] The fre- 
quency range considered is from nearly zero up to slightly 
less than twice the Nyquist frequency. The maximum Lomb 
power is 6.10 and it occurs at a frequency of 8.47 cycles/year. 

A simple way to assess the significance of a peak in such 
a spectrum is to make a histogram of the number of frequen- 
cies as a function of power. In the absence of any time varia- 




0.0 2.0 4.0 6.0 



Power 

FIG. 6: Histogram of powers in spectrum of Fig. |5] The bin size is 
0.07 power units. The solid line is the expected distribution if there 
is no time variation, i.e., the number of frequencies x exp(-power), 
integrated over the limits of each bin. 



tion this distribution is an exponential; if there were any peak 
present with significant power it would appear at the upper 
end of the distribution and be clearly separated from the ex- 
ponential trend. This distribution for the spectrum of Fig.|5]is 
shown in Fig. |6] As this distribution visually shows, there is 
no evidence for exceptionally high power in the data spectrum 
at any frequency. 

A quantitative way to prove that no frequency has excep- 
tionally high power is with a shuffle test. In this test the 
SNU results are randomly re-assigned to the different runs, 
the power spectrum is recalculated, and the maximum power 
is found. The maximum power in the spectra from 1000 such 
shuffles is plotted in Fig. [7] The observed maximum power 
for the SAGE data of 6.10 occurs very near the center of this 
distribution. Of the shuffles, 44% have a greater power than 
for the data and 56% a lesser power, thus showing that the ob- 
served power distribution is consistent with the assumption of 
a constant rate. 

The same result is obtained if power spectra are produced 
over a wider frequency range. For ranges up to 19 per year and 
38 per year the maximum power remains 6.10 at 8.47 cycles 
per year and is not statistically significant. 

The three tests in this section do not find any evidence for 
peiiodic variation in the SAGE data. The frequency range to 
which these tests are sensitive is about one cycle per month 
to about one cycle per ten years. For a frequency to be de- 
tected within this range the amplitude of the peiiodic oscil- 
lation must, as shown by Pandola |37|, exceed the statistical 
uncertainty of a single run, i.e., must be > 30 SNU. 

We end this section by noting that the winter minus summer 
difference in SAGE capture rate is R^f^ - Rs - 5.8+gj SNU 
where the stated unceitainty is only statistical. For this calcu- 
lation summer was defined as the +l4-year interval centered 
on 2 1 June and winter as the rest of the year. In our method of 
data analysis f2] we remove the known change in rate caused 
by the Faith's orbital eccentricity. If, rather than using the 
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FIG. 7: Histogram of maximum power found in Lomb power spec- 
trum analysis of 1000 random shuffles of the 168 SAGE runs. The 
bin size is 0.13 power units. 

above solstice-based definition, we define summer as the ±14- 
year interval centered on 5 July, the time of the aphelion, then 
R^-Rs = 4.2+^-2 SNU. With both of these definitions /Jw-^s 
is consistent with zero, indicating that there is no appreciable 
difference between the day and night capture rates in Ga, as is 
expected for the currently determined values of the neutrino 
oscillation parameters 13S\ ). 



VIII. SUMMARY AND DISCUSSION 

In 18 years of operation the SAGE experiment has car- 
ried out 168 measurements of the solar neutrino capture rate. 
Analysis of the times of occurrence of all events in windows 
centered on the L and K peaks and with 2 FWHM width as- 
cribes 853.9 of these events to the decay of ^'Ge. To com- 
pare these values with other radiochemical solar neutrino ex- 
periments, the CI experiment collected solar data for a longer 
time (108 runs with rise-time counting during 24 years), but 
the number of detected events within an energy window of 
2 FWHM that were ascribed to ^'^Ar was 875 fl |, quite com- 
parable to what we now have in SAGE. So far as we are aware, 
neither Gallex nor GNO have reported their total number of 
detected ^'Ge events, but because the mass of their Ga target 
was less (30 tonnes), and their number of runs was less (123 
during 12 years), this number must be considerably less than 
we now have in SAGE. 

The measured best-fit capture rate in the SAGE experiment 
is 65.4+3Q (stat) (syst) SNU. Combining this with the re- 
sults of Gallex and GNO gives the rate as 66.1 + 3.1 SNU, 
where statistical and systematic uncertainties have been com- 
bined in quadrature. The solar-model prediction for the Ga 
experiment is 66.3 SNU for a model with high metallicity and 
63.2 SNU for a model with low metallicity, where the un- 
certainty of both estimates is ~4%. There is thus excellent 
agreement between theory and experiment for the Ga exper- 
iment. Further, both the experimental measurement and the 



theoretical prediction are known to about the same accuracy. 

By use of the results of other solar neutrino experiments 
and neutrino oscillation theory we have derived the contem- 
porary value of the pp flux from the Sun to be (3.40^^4^) x 
10"V(cm2 s) at the Earth and (6.0 + 0.8) x lO'^Acm^'s) at 
the Sun. The latter is in good agreement with standard so- 
lar model predictions of 5.97 ± 0.04 (high metallicity) and 
6.04 + 0.03 (low metaUicity), both in units of lO"* vj{cm^ s). 
Gallium experiments have thus proven that the overwhelming 
fraction of solar neutrinos that reach the Earth are the low- 
energy neutrinos from the pp reaction. 

We have assumed in these calculations that the cross sec- 
tion for capture to the two lowest-lying excited states in ^'Ge 
is zero, as is implied by the four neutrino source experiments 
with gallium. This assumption is in contradiction to the stan- 
dard interpretation of the two experiments that have attempted 
to measure the Gamow-Teller strength of these low-lying ex- 
cited states in ^'Ge. These experiments were made by (p,n) 
scattering |Il9l|39l and (^He,f) scattering |40 41 , 42 1. If all the 
events observed in these experiments at low excitation energy 
are attributed to Gamow-Teller strength the results of these 
experiments are in reasonably good agreement. It is, however, 
not evident that these experiments solely measure Gamow- 
Teller strength-as emphasized by Haxton |,18J , for very weak 
transitions, such as is the case in these experiments, there may 
be an appreciable (perhaps dominant) contribution to the cross 
section from the spin-tensor interaction. New experimental 
data is needed to settle this question and to definitively deter- 
mine the magnitude of the matrix elements for neutrino cap- 
ture to these two low-lying excited states. We strongly en- 
courage any new experiments that might shed light on this 
question. As part of our future experimental program we in- 
tend to pursue a new measurement that will use a very intense 
neutrino source in an optimized detector geometry. 

The SAGE experiment continues to collect data on the so- 
lar neutrino capture rate with a gallium target. Up to now it 
is only the Ga experiment that has measured the low-energy 
pp solar neutrinos. As we continue to monitor the solar neu- 
trino flux we will increase our statistical accuracy and further 
reduce our systematic uncertainties. 
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APPENDIX A: CALCULATION OF EXTRACTION 
EFFICIENCY FROM ISOTOPIC ANALYSIS 

It is assumed the extracted Ge consists of a combination of 

• Ge from the carrier added for the current extraction, 

• residual Ge that remained from the carrier added for the 
two preceding extractions, and 

• additional Ge with natural isotopic composition, such 
as may be dissolved from the surfaces of the extraction 
system or the vessels that contain the Ga. 

According to this model the predicted mass of isotope ; ob- 
tained in extraction «, called Mj^(i), is thus 

M^(0 = e„{C„(0+(l-e„-i)[C„_,(0+(l-e„-2)C„-2(0])+£„/(0, 

(Al) 

where e,„ is the efficiency of Ge removal in extraction m, Cm{i) 
is the mass of isotope / of carrier added to extraction m (where 
m can take on the values n,n—\,orn- 2), and £„ is the mass 
of additional Ge with natural isotopic composition that is 
removed in extraction n. It is assumed that /(/), C„{i), Cn-\{i), 
and Cn-iii) are known and the variables e„, e„_i, e„_2, and E„ 



are to be determined. For extractions there are thus 2N 
variables (the extraction efficiency and mass of extra natu- 
ral Ge for each extraction) and 5A^ equations that relate these 
variables, one for each of the naturally-occurring Ge isotopes, 
^"Ge, ^^Ge, ^^Ge, ^"^Ge, and ^^Ge. Because there are more re- 
lationships than unknowns, the problem is solved by finding 
the set of variables that minimizes the function 



^2 ^ ^ ^l MP{i)-Ml{i) 



(A2) 



where M^(i) is the measured mass of isotope / in extraction 
n and cr„(/) is the total uncertainty in the knowledge of the 
predicted and measured masses in extraction n. 



APPENDIX B: RESULTS FOR EACH SAGE EXTRACTION 

The capture rate for each SAGE extraction is given in Ta- 
ble |v] The statistic Nw^ in this table measures the goodness 
of fit between the observed sequence of events in time and the 
time distribution predicted by the model used in analysis, viz., 
that the events are produced by the sum of two processes: the 
decay of a fixed initial number of ^'Ge atoms and background 
events at a constant rate. The probability to obtain a value of 
Nw^ larger than in the observed time distribution is given in 
the last column. It is derived by 1000 simulations for each ex- 
traction using the method in Ref. |35 | and has an uncertainty 
of ~l.5%. The time that the exposure interval for each ex- 
traction began and ended can be obtained from the entries in 
columns 2 and 3 by use of transformation equations (6.4) in 
Ref. 121. 



TABLE V: Results of analysis of SAGE extractions. 
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APPENDIX C: MODIFIED CROSS SECTION 
FOR NEUTRINO CAPTURE 



The cross section for neutrino capture by ^'Ga was calcu- 
lated by Bahcall |19| and is given in his Table II (best es- 
timate). Table III {3cr lower limit), and Table IV (3cr upper 
limit). Based on information on the contributions of the ex- 
cited states given in the text of Bahcall's article, if we assume 
the matrix element for neutrino capture to the first two excited 
states of ^'Ge to be zero, but that the matrix elements of the 
other excited states are unchanged, we can approximate the 
best-estimate cross section in various energy regions as fol- 



lows: 



CT = 



13.10 H-91.29(£v 

0.946 CTbest 
0.953 CTbest 
0.96 CTbest 
0.99 (Tbest 



.Q 24)1.157 for 0.24 <£^< 0.733 
for 0.733 < 1.033 
for 1.033 < 1.483 
for 1.483 < £,< 1.983 
for 1.983 < £v< 30.0 



where cTbest is the value in Table II of Ref. (19], cr is in 
10 cm^, and the neutrino energy Ey is in MeV. The re- 



sults are given in our Table VI The ±lcr limits in Table VI 
were obtained in a similar manner from the 3cr limits given in 
Ref |,19J. 
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TABLE VI: Approximate cross section for neutrino capture by ^'Ga if the contribution of the first two excited states is set to zero. 





Cross section (10 


cm") 


V energy 


Cross section (10 


cm^) 


V energy 


Cross section (10 **' 


cm^) 


(MeV) 


Best 


-la- 


+ la" 


(MeV) 


Best 


-la- 


+ ltT 


(MeV) 


Best 


-la- 


+ 1(T 


0.240 


1.310X 10' 


1.280 X 10' 


1.340 X 10' 


1.445 


1.944 X 10^ 


1.897 X 10^ 


2.274 X 10^ 


9.500 


4.749 X 10* 


4.053 X 10* 


6.275 X 10* 


0.250 


1.357 X 10' 


1.326 X 10' 


1.388 X 10' 


1.500 


2.153 X 10^ 


2.104 X 10^ 


2.495 X 10^ 


10.000 


5.653 X 10* 


4.820 X 10* 


7.474 X 10* 


0.275 


1.499 X 10' 


1.465x10' 


1.533 X 10' 


1.600 


2.451x10^ 


2.395 X 10^ 


2.859 X 10^ 


10.500 


6.638 X 10* 


5.650 X 10* 


8.785 X 10* 


0.300 


1.662 X 10' 


1.624 X 10' 


1.700 X 10' 


1 .700 


2.771 X 10^ 


2.707 X 10^ 


3.252 X 10^ 


1 1 .000 


7.703 X 10* 


6.548 X 10* 


1.020 X 10= 


0.325 


1.836 X 10' 


1.794 X 10' 


1.878 X 10' 


1.750 


2.939 X 10^ 


2.871 X 10^ 


3.458 X 10^ 


1 1.500 


8.848 X 10* 


7.510 X 10* 


1.173 X 10= 


0.350 


2.018 X 10' 


1.972 X 10' 


2.064 X 10' 
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3.932 X 10^ 
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4.712 X 10^ 
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1.007 X 10^ 


8.535 X 10* 
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6.428 X 10^ 


5.986 X 10^ 


7.826 X 10^ 
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